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S~ escaping and scattering probabilities in the (K~ ,K + ) reaction with nuclear targets 
is investigated on the basis of the eikonal approximation. Calculations are performed using 
the S = —2 sector of the Nijmegen model-D interaction, whose hard core radii are treated 
phenomenologically. The S~ N elastic cross sections are derived from the 9 He(K~ , K + ) data 
of the BNL-E906 experiment, which is useful to determine the hard-core radii in P-states. 
The S single-particle potential energy is investigated with the G-matrix calculation in nuclear 
matter, which is sensitive to the hard-core radii in S-states. The important constraint for 
.— i. , the EN interaction is obtained by combining the two analyses. 

r-> 

§1. Introduction 

It is very important to study the properties of SN and S'-nucleus interactions 
in hypernuclear phenomena, because free-space SN scattering experiments are dif- 
ficult to perform with current experimental apparatus. Interesting information was 
obtained from events of simultaneous^missions of two A hypernuclei (twin A hyper- 
nuclei) in the KEK-E176 experiment. GP The obtained values of binding energies B s - 
lead to an important conjecture for the S'-nucleus potential. The BNL-E885 collab- 
orations measured the missing mass spectra for the 12 C(K~, K + )X. reaction. Rea- 
sonable agreement between this data and theory is realized by assuming a ^-nucleus 
potential U~(p) = —Vof(r) with well depth Vq ~ 14 MeV within the Wood-Saxon 
(WS) prescription. 

Recently, a new experiment (BNL-E906) has been performed in order to detect 
charged particles emitted from (K~ , K + ) reaction points with a 9 Be target using a 
cylindrical detector system. eP Our concern in this work is with the emission probabil- 
ities of S~ particles produced through the quasi-free process in the target. In recent 
analysis,© S~ emission events have been separated into two categories, the S~ es- 
cape process without any interaction after its production, and the scattering process 
with a target nucleon, for which the obtained probabilities are 63.6 ± 8.2% (escap- 
ing probability) and 14.5 ± 2.6% (scattering probability), respectively. Prom the 
scattering probability, the S~N elastic scattering cross section in a nuclear medium 
was found to be 20.9 db 4.5(stat)l 2 ' 4 (syst) mb by using the eikonal approximation, 
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in which the number of the nucleons was eight and the reaction cross sections of K~ 
and K + with the nucleons were 28.0 mb and 19.5 mb, respectively. Also, the elastic 
scattering probability was separated into the E~p part (5.8 ± 3.5%) and the E~n 
part (8.7 ± 4.4%). 

Nara et al. studied the K + momentum spectrum from (K~,K + ) reactions in 
the framework of inter-nuclear cascade-type calculations.^ Here, two-step processes 
such as K~p — > M'Y followed by M'N — > K + Y, where M' represents intermediate 
mesons, play a significant role in K + yields, especially for heavy targets. The contri- 
butions from these processes seem to be non-negligible, even in the high-momentum 
region of K + , where quasi- free E~ productions are dominant. Because the E~ emis- 
sion probabilities are defined as the ratio of the number of emitted E~ particles to 
produced E~ particles, the two-step contributions to (K~,K + ) events make them 
larger. The ratio of the two-step processes in the (K~ , K + ) events, denoted as a tW oj 
is estimated to be about 10% at most.Er The above values of the escaping proba- 
bility, 63.6%, and the scattering probability, 14.5%, for a two = 0% change to 67.1% 
(71.0%) and 15.1% (15.8%), respectively, for a two = 5% (a two = 10%). 

In the report of the E906 experimentH^ the most important results of the theo- 
retical calculations are cited, LP'EP'EiP where a simple model is proposed for the E~ 
emission process. The outline of this paper is as follows. First, the plausibility of this 
model is investigated in light of the E906 data. Second, the observed E~ escaping 
and scattering probabilities are shown to give an important constraint for the EN 
interaction, which complements the E well depths derived from the data of E885 
experiment. 

§2. E~ emission probabilities and the E ~N elastic cross section 

We now formulate our model. LP'tP'tO) The transition rates w s £ at due to elastic 
scattering between a quasi- free E (momentum p=) and nucleons (momenta pjv) in 
nuclear medium are defined by 

wT = J d 3 p N J d 3 ff SN J d 3 P' 3N p f N (p N ) v SN ^ 

x S((P~N ~ p'sn)/p'In) * 3 (Pzn - P'sn) Q(Pn - PfWb - qs), (2-1) 

where p is the nuclear density, 4j| is the EN differential scattering cross section for 
the transition from the initial relative momentum pen to the final one p'^jy, and v~n 
is the relative velocity. Also, Pen and P^n denote the initial and final center-of-mass 

momenta, respectively. The function /n(pn) = (I^Pf) 6{pf — Pn) represents a 
Fermi distribution of nucleons, with pp being the Fermi-momentum. The factor 
9{p' N ~pf) is the Pauli exclusion factor for the final nucleon state of momentum p' N , 
and 0(p' s — qz) allows transitions only to E states whose momenta (p~) are larger 

than the critical momentum qs determined by ^jffz + U~(p) = 0, U~{p) being a E 
single-particle potential in medium. This constraint for E states results from the 
fact that lower-momentum components of scattered E particles are considered as 
being absorbed into target nuclei and not detected experimentally. In this work we 
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assume Ug(p) = Us{po){p/ Po)-, where po is the normal density. The value oiJJs(po) 
is given roughly by the Wood-Saxon well depth of E in the E885 experiment. EP Here, 
we use Us(po) = — 14 MeV, but our results are not sensitive to the value of U~(po)- 

Assuming isotropic angular distributions tjj = ~£f") Eq. (2-1) can be reduced to a 
compact form, OP where ct^ n is the angle-integrated SN elastic cross section. In the 
above expressions, charges of S particles and nucleons are implicit for simplicity. In 
actual calculations, we treat explicitly the processes S~p — ► S~p and S~n — ► S~n. 
Then, the transition rates of S~ for protons and neutrons in the target are denoted 
as w^^ip) and w s ^(n), and we have w B £^ = w^^(p) + w B £^(ri). 

We define here the E~ absorption rate w^O. phenomenologically This is deter- 
mined from the probability of S" vanishing due to processes such as S~p — » S°n 
and S~p — > AA in target nuclei. The corresponding S~ absorption cross section 
is denoted as a a bs- The S~ mean free path (MFP) Act- in medium is given as 
follows. The MFPs for scattering and absorption cross sections are represented by 
A?S- = v s - /w s ^ and A*.! = l/(/90" a bs)j where uc— is the velocity of S~ . Then we 
have the relation 1/Act- = l/Afi- + 1/A^f . Using Eq. (2-1), Af- is obtained as a 
function of pg- and p. 

The probability of E~ reaction processes in a finite nucleus can be estimated on 
the basis of the eikonal approximation. Assuming forward scattering, 9^+ = 9s- = 
0°, the (K~,K + ) effective number for the reaction (scattering and absorption) of 
the produced quasi-free S~ and nucleons is given by 



Af^ c =/ 27rb db / dzp(Vb 2 + z 2 ) 

JO J-oo 



>0 J-oo 

x exp< —a K - N / p(Vb 2 + z' 2 )dz' - a K + N 



oo 




x r exp (-f^')r 

with 1/Act- = I/A??-* + 1/A^ 3 - • Here, cr K - N and o~ K + N are the isospin averaged 
total cross sections between if - and K + and a nucleon. We set cr K - N = 28.0 
mb and o~k+n = 19.5 mb in this work. The total effective number A tota i for the 

(K~, K + ) reaction is given by omitting the factor si — expf — J^° j — dz' I > in the 

above expression, and the reaction probability of E~ in medium is given by 

prea, = Ngc/N^ . ( 2 . 3 ) 

Similarly, we obtain the effective number JV-„ for Z~ absorption in medium by 
replacing Act- in Eq. (2-2) with A*-- Then, the scattering part of the effective 
number, which originates from the S~N scattering in medium, is given as 

/■OO /"OO 

NW- t = 2irbdb dzp{Vb 2 + z 2 ) 

JO J-oo 

x exp< —a K - N / p(Vb 2 + z' 2 )dz - a K + N / p( \Jb 2 + z' 2 )dz' > 
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x ^-f ^{'-^K-f at* 2 ')} • (24) 

Then the E~ scattering probability for the E~ emission after E~N scattering in 
medium is given by 

Pf at = Ar-at/^^^ _ (2 . 5) 

The E~ escaping probability for E~ emission with no interaction after the production 
is given by 

Pf- C = 1 - A^/AWl • (2-6) 

The E~ emission probability is given as the sum of the above two probabilities: 

pcinit -i pabs pesc _i_ pscat ('}.7\ 

with the E~ absorption probability PJ. 1 ! 8 = N^l /-/Vtotal- The probabilities PJPf, 
pscat ^ pemit an( ^ pabs can ^ e CO mpared with the E906 data. Here, it is very impor- 
tant in our treatment that the experimental values of P~f and PSFf 1 * be obtained 
separately in the analysis of the E906 data. As another expression for the E~ scat- 
tering probability, we define here the in-medium scattering cross section cr S cat by the 
relation A^-* = ^ / (po'scut) ■ Here, the value of a sca t can be obtained phenomenologi- 
cally without using Eq. (2-1), by treating it as a parameter in Eq. (2-4) to reproduce 
the data of Pf^. 

The proton and neutron distributions in 9 Be are represented as p(r) = po(l + 
a(r/1.65) 2 ) exp(— (r/1.65) 2 Lwith a=1.46 (proton) and 2.10 (neutron) on the basis 
of a theoretical calculation. eP 

It is necessary to use appropreate EN interactions for calculations, though at 
present we do not have sufficient knowledge of realistic models. In this work we use 
mainly the S = —2 sectors of S*C/(3)-invariant OBE models, which have been used 
successfully to explore properties of S = — 2 hypernuclear systems such as double- 
A nuclei and S'-micleus bound states. Among the various interaction models, the 
Nijmegen D (ND)liiP seems to be suitable at present, because this model reproduce 
the attractive S" _ -nucleus potentials found in the KEK-E885 experiment a> with a 
reasonable choice of the hard-core radii r c , as shown below. ND should be considered 
here as a useful model for our purpose to analyze the data of both the S~N cross 
section and the E well depth in medium complementarity. Another ■S'./V" interaction, 
which we use only for comparing with ND, is of a purely phenomenological Wigner- 
type represented by a one-range Gaussian (ORG) form v{r) = uoRGexp(— (r//3) 2 ) 
with (3 = 1.0 fm, for which the isospin dependence is not taken into account. 

Now our concern is in determining how the E906 data can be used to determine 
hard-core radii r c in ND (t>oRG i n ORG). Hereafter hard-core radii in S-, P- and in- 
states are denoted r c (S), r c (P) and r c (D). The E906 data gives the EN scattering 
in medium at E momentum around 550 MeV/c, where the P-state contributions are 
more important than the 5-state contributions, and the results are sensitive to the 
value of r c (P), not only to that of r c (S). Contrastingly, the S-state contributions are 



3 emission probabilities at (K ,K + ) reaction points 5 

Table I. In the cases of awo = 0.0, 0.05 and 0.10, the values of cr a b s (nib) and r c (P) (fm) are chosen 
so as to reproduce the experimental values of P^I, Pgf and PJPf*. The notation R(2~p/E~n) 
is used to represent P| c f t (p)/P|Pf t (n). The values of PJ!f are Pgpf* are the fitted ones and agree 
with the experimental values. See the main text for the definitions of cr a b s (mb) and <7 scat (mb). 
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far more important than the P-state contributions for 3 well depths in nuclei. In the 
case of ND, therefore, the two data can be used, so that the hard-core radii of S- and 
P-states are determined separately. First, the value of r c (S) = 0.50 fm is chosen so 
as to result in an appropriate attraction of the calculated 3 single-particle potentials, 
as shown in the following section. Next, the value of r c (P) is determined so as to 
reproduce the E906 data reasonably. Additionally, we assume r c {D) = r c {S) = 0.50 
fm for simplicity. The isospin dependence of the hard-core radii of ND is not taken 
into account. It should be noted, however, that the isospin dependence in the OBE 
part of ND is treated precisely. In the case that we use ND, the 3N — A A coupling 
interaction is treated exactly, but the couplings to AS and EE channels are ignored 
for simplicity. 

We now give our results. There are essentially two free parameters in our formal- 
ism, o" a bs an d r c (P) in ND (^org m ORG). The former is chosen so as to reproduce 
the experimental value of P2r s , namely the probability of vanishing of 3~ in the 
(K~,K + ) reaction. We have P| b _ s = 0.219, 0.178, 0.133 O abs = 21.1 mb, 16.4 mb, 
11.7 mb) for a two = 0.0, 0.05, 0.10, respectively. Next r c (P) in ND (v RG in ORG) 
is determined so as to reproduce the experimental values of P~£ and P-i 1 *. Table I 
gives the calculated values with ND at £>=■- = 550 MeV/c, around which the peak 
value of the 3~ production strength is obtained. Here, the experimental values of 
para ^ pesc an( j pscat ^ urn OUTj to De reproduced precisely with the values of <7 a b s and 
r c {P) given in the table. The in-medium scattering cross sections <7 scat are deter- 
mined here so as to reproduce the data of P2F_ phenomenologically by using the 
relation A~_ = l/(/0<7 scat ), instead of deriving A~_ from the 3N interactions. 

In the E906 experiment, the elastic scatterings between the produced 3~ par- 
ticles and nucleons in the target were separated into the 3~p part (5.8 ± 3.5%) and 
the 3~n part (8.7 ± 4.4%). d' The ratio of the probabilities of 3~p and 3~n scat- 
terings are regarded as corresponding to R(3~p/3~n) = Pl cat (p)/P2, cat (n) in our 
model. The calculated values are given in Table I. These are found to be similar 
to the experimental value, ~ 0.67. It is more reasonable, of course, to also consider 
the isospin dependence of r c (S) simultaneously. There is no information, however, 
allowing us to determine them independently at present. 

Similar results are obtained with ORG. For instance, the above values of P2?^ and 
pscat £ or at ^ o = o.o can be reproduced with i>org = —54.5 MeV. Because no isospin 
dependence is taken into account in ORG, the obtained value of R(3~p/3~n) = 
0.603 is significantly different from that with ND. 

In Fig. 1, the calculated values of the 3~ escaping probabilities P2?- and the 
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Fig. 1. Calculated values of P* B f and Pj^** in the case a two =0 as functions of p s - in (a) and (b). 
Here, the solid and dashed curves correspond to ND and ORG. The experimental values of 
Pff = 63.6 ± 8.2% and P| rat = 14.5 ± 2.6% at p s - = 550 MeV/c are shown for comparison, 
where the dotted horizontal bars indicate the width of the S~ spectrum (FWHM). 



E scattering probabilities P*Ff* are plotted as functions of pg- in (a) and (b), 
respectively. Here, the solid and dashed curves correspond to ND and ORG. The 
experimental values of Pf s f = 63.6 ± 8.2% and P| cat = 14.5 ± 2.6% at p E - = 550 
MeV/c are shown for comparison, where the dotted horizontal bars indicate the width 
of the ~- spectrum (FWHM). The values of Pff and Pf at for ND and ORG turn 
out to be adjusted so as to reproduce the experimental values at pg- = 550 MeV/c, 
though their pg- dependences differ significantly. Thus, it is found that our two EN 
interactions result in fits of equal quality by adjusting their interaction parameters. 
As demonstrated in the next section, we need to complement arily analyze the other 
data, such as those for the E well depths, in order to make clear the features of the 
respective EN interactions. 

In Fig. 2, the calculated values of the E~N elastic cross sections are plotted 
as functions of the E~ momentum pg- (MeV/c) for ND with r c {P) = 0.480 fm, 
0.491 fm and 0.502 fm, by solid, dashed and dotted curves, respectively. The E~N 
cross sections were obtained by averaging the E~p and E~n cross sections with the 
weight of 3 protons and 5 neutrons. The ambiguity concerning the two-step processes 
is not serious for the resulting E~N cross sections. This is essentially because E~N 
scattering events were observed explicitly in the E906 experiment, in spite of the 
unknown two-step contributions. 

There is an experimental error of 2.6% for PSf* = 14.5% in the case that we 
assume a t wo=0; explicitly, the upper and lower values of P~£ are 17.1% and 11.9%. 
Here, we attempt to reproduce the deviation of 2.6% by changing only the value of 
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Fig. 2. Calculated values of S N elastic cross sections in mb are plotted as functions of the 5 
momentum p s - (MeV/c) for ND with r c (P) = 0.480 fin, 0.491 fm and 0.502 fm, by solid, 
dashed and dotted curves, respectively. The S~ N cross sections are obtained by averaging the 
E~p and E~n cross sections with the weight of 3 protons and 5 neutrons. 



r c (P): The values of 17.1% and 11.9% are reproduced by choosing r c (P) = 0.525 
and 0.447, respectively. This choice is expressed as r c (P) = 0.480^ ; 03 3 fm for 
convenience. In the above analysis, we assumed that the values of at wo is between 
and 10%, giving rise to r c (P) = 0.480 - 0.502 fm. Thus, it turns out that the 
ambiguity due to the two-step process is less than the experimental error bar. 

§3. 3N G-matrix calculation and Us 



In order to demonstrate the properties of our interaction, let us perform the 
3N-AA channel coupled G-matrix calculations in symmetric nuclear matter with 
ND whose hard-core radii are taken as the above values. Calculations here are done 
in the framework used in a previous workEj-*: We adopt the simple QTQ prescription 
for the intermediate-state spectrum, which means that no potential term is taken 
into account in the off-shell propagation. 

In the previous analysis with ND, r c (S) and r c {D) are fixed at 0.5 fm, and r c (P) 
is adjusted so as to reproduce the E906 data. Because our determined values of r c (P) 
are around 0.5 fm, it is reasonable here to use the simple choice r c (S) = r c (P) = 
r c (D) = 0.5 fm. In this case, we calculate the single-particle potential energies Us for 
a zero-momentum 3 at normal density. The result is Ug = —14.3 MeV, in which the 
5-state, P-state and .D-state contributions are Us(S) = 2.8 MeV, Us(P) = —15.6 
MeV and Us(D) = — 1.6 MeV, respectively. It should be noted here that the negative 
value of Us is due to the strongly attractive P-state contribution. As shown in the 
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previous section, P^ = 14.5±2.6% corresponds to r c (P) = Q.A&Q'^'^m, which leads 
to Us = — 15.1T 1 ' 8 MeV. Thus, the ambiguity of Us due to the experimental error 
of the E906 data is found to be small. 

We now demonstrate that our obtainedvalues of Us are consistent with the ex- 
perimental results. In the E885 experiment, o* the missing mass spectra for l2 C(K~ , K + )X 
reaction was measured. Reasonable agreement between the data and theory is 
achieved by using a Wood-Saxon (WS) potential between S~ and the n B core 
whose well depth is about 14 MeV. In order to compare our result with this WS 
potential, we construct the S- nucleus potential from the calculated values of Us(p) 
with a local density approximation (LDA) as 



U s ( 



(t0F)~ 3 [d7U s (p(r'))exp(-\r-7\ 2 /t 2 )p(r') , (3-1) 



where the_finite-range effect of SN interactions is taken into account by the pa- 
rameter tX-3) The nuclear density distributions p(r) are obtained from the Skyrme- 
Hartree-Fock calculations. We compare our obtained S'- 11 B potential with the above 
WS potential by calculating the volume integral per nucleon, Jy /A = — A^ 1 J U(r)d 3 r, 

the root mean square radius, (Ry) 1 ' 2 = [J U(r)r 2 d 3 r/ J U(r)d 3 r] , and the bind- 
ing energy of the ground IS" state, -B~(15), obtained without taking the Coulomb 
interaction into account. 

In the case of the above WS potential, the obtained values are Jy j A = 133 MeV 
fm 3 , (R v } 1/2 = 3.1 fm and £-(15) = 2.2 MeV. By setting t = 1.0 fm in the case 
of our G-matrix interaction derived from ND with r c (5) = r c {P) = r c (D) = 0.5 fm, 
we obtain Jy j A = 126 MeV fm 3 , {R v ) l/2 = 2.9 fm, and 5= (15) = 2.2 MeV. Thus, 
it is confirmed that the two ^-nucleus potentials are very similar. This results from 
our choice of r c (S) = 0.5 fm, because changes of r c {P) and r c (D) have far weaker 
effects on the Us values. 

The important information regarding the S'-nucleus potential can be obtained 
from the events of simultaneous emissions of two A hypernuclei (twin A hypernuclei) 
in the KEK-E176 experiment. These data give the energy difference between the 
initial S~ bound state and the final twin A state, namely the binding energy -B=— 
between S~ and the nucleus. For instance, the event studied in Ref. |l|) indicates the 
2P ~~ bound state in 12 C with B S (2P) = 0.57 ± 0.19 MeV. With the above proce- 
dure we can derive the S~- 12 C potential straightforwardly. Including the Coulomb 
interaction, we obtain the calculated value of 0.40 MeV for Bs{2P). This is slightly 
smaller than the experimental value, but it is still within the error bar. 

Finally, let us point out the important contributions of the ^A^ P-state inter- 
action to the single-particle potential energy of S in dense nuclear matter, which is 
an interesting feature of ND. Considering the important roles of negatively-charged 
baryons in neutron stars, we calculate here the s.p. potential Us- in neutron mat- 
ter. EJ) In Fig. 3, the Us- values for ND with r c (S) = r c (P) = r c (D) = 0.5 fm are 
plotted as a function of p/po by the solid curve, where po = 0.17 fm~ 3 is the normal 
density. The S-, P- and D-state contributions are plotted by the dashed, dotted 
and dot-dashed curves, respectively. It should be noted here that the 5-state contri- 
bution is strongly repulsive and the negative values of Us- are due to the strongly 
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Fig. 3. The values of U B - in neutron matter calculated with ND with r c (S) = r c (P) = r c (D) = 0.5 
fm are plotted as a function of p/po by the solid curve. S-, P- and D-state contributions are 
plotted by the dashed, dotted and dot-dashed curves, respectively. 



attractive contribution of the P-state interaction. In the Ref. 15), similar calcula- 
tions were performed using ND with hard-core radii assumed to have the same values 
as the corresponding ones in the NN channel. The value of r c (P) = 0.346 fm in 
that case is significantly smaller than the value r c (P) ~ 0.5 fm in the present case. 
This implies that the former P-state interaction is too attractive. 

§4. Conclusion 

The BNL-E906 experiment was performed in order to detect charged particles 
emitted from (K~ , K + ) reaction points on 9 Be target. In recent analysis of the S~ 
emission events, the S~ escape process without any interaction after its production 
and the scattering process with a target nucleon have been extracted separately. 
These are specified by escaping and scatteringjprobabilities, respectively. A simple 
model was proposed for these processes, LP'EP'tJ) which makes it possible to relate 
these probabilities with underlying EN interactions. 

Calculations were performed using the S = —2 sector of ND, where the hard- 
core radii are treated as free parameters. The obtained S~ escaping and scattering 
probabilities are consistent with the E906 data, where the choice of hard-core radii 
in the P-states is sensitive to the results. 

Another important piece of information is obtained from the E single-particle 
potential energy, which is derived from the EN G-matrix calculations with ND in 
nuclear matter. The resulting value is very sensitive to the hard-core radii in the 
5-states. Our obtained S'-nucleus potential, constructed with LDA, is similar to the 
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WS potential obtained from the E885 experiment, 

Thus, data from the E906 and E885 experiments can be used complement arily, so 
that the hard-core radii of ND in S- and P-states can be determined separately. The 
determined values of hard-core radii are not so different from those in S = and S = 
— 1 channels. This fact indicates that ND is a useful interaction model also in S = —2 
channels, and we expectthat it will be helpful to explore the effect of S~ mixing 
in neutron star matter. i3 Of course, there still remain some problems regarding 
modeling of ND. Similar calculations are needed with use of more sophisticated 
interaction models. The present study given will be an important test for such EN 
interaction models. 
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